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ABSTRACT 

Turbine vane heat transfer predictions are given for 
smooth and rough vanes where the experimental data 
show transition moving forward on the vane as the sur- 
face roughness physical height increases. Consistent with 
smooth vane heat transfer, the transition moves forward 
for a fixed roughness height as the Reynolds number in- 
creases. Comparisons are presented with published exper- 
imental data. Some of the data are for a regular rough- 
ness geometry with a range of roughness heights, Reynolds 
numbers, and inlet turbulence intensities. The approach 
taken in this analysis is to treat the roughness in a statis- 
tical sense, consistent with what would be obtained from 
blades measured after exposure to actual engine environ- 
ments. An approach is given to determine the equivalent 
sand grain roughness from the statistics of the regular ge- 
ometry. This approach is guided by the experimental data. 
A roughness transition criterion is developed, and com- 
parisons are made with experimental data over the entire 
range of experimental test conditions. Additional com- 
parisons are made with experimental heat transfer data, 
where the roughness geometries are both regular as well 
as statistical. Using the developed analysis, heat transfer 
calculations are presented for the second stage vane of a 
high pressure turbine at hypothetical engine conditions. 

N omenclat ur e 

A + Near wall damping coefficient 

C Modeling constant 

(7x Axial chord 

Cf/2 Friction coefficient 

d Base diameter 

H+ Normalized roughness height 

h Heat transfer coefficient 

K Acceleration parameter 

k Physical roughness height 

ks Equivalent sand grain roughness 

M Mach number 

n Number of points in roughness calculation 

N Spot production parameter 

Nu Nusselt number 

P Pressure 


p Pitch 

R & Absolute roughness height 

Rrms Root mean square roughness height 
Re Reynolds number 

Re i Unit Reynolds number 

5 Surface distance 

Sk Skewness 

St Stanton number 

Tu Turbulence intensity 

U Velocity 

Y Loss coefficient 

y Measured height 

Distance from surface 
y Average height 

Ay Increase in mixing length due to roughness 

Ay + Normalized value for Ay 

6 Momentum thickness 

n Mixing length constant 

p Dynamic viscosity 

/i t ,i Inner region eddy viscosity 

p Density 


Subscripts 

2 

FS 

IN 

RLM 

ROUGH 

SMOOTH 

SK 

SL 

ST 

t 


Vane exit 
Local freestream 
Vane inlet 
Relaminarization 
Rough 
Smooth 
Skewness 
Length 

Transition start 
Total 


INTRODUCTION 

Surface roughness adversely affects turbomachinery 
performance by increasing external heat transfer, and by 
increasing loss. Measured surface roughness for in service 
blades show a high degree of variability. It is useful to 
the designer to have an estimate of the effects of surface 
roughness on both heat transfer and aerodynamic perfor- 
mance. Data show that as the roughness height progres- 
sively increases, turbine vane transition moves forward on 
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both the suction and pressure surfaces of the vane. Data, 
such as those of Arts[l], show that a smooth vane can have 
laminar flow over much of the vane at moderate Reynolds 
numbers, even with relatively high inlet turbulence. When 
the boundary layer is not forced to be turbulent due to film 
cooling, understanding the interaction of surface roughness 
and transition is important. The change in blade surface 
heat transfer with transition is a very good indicator of 
transition start and length. Tests of a solid high pres- 
sure turbine vane are therefore valuable, even if the vane 
is film cooled when installed in an engine. Transition is 
typically not an issue for the first turbine stage, since the 
blades are generally film cooled. Arts [2] has shown that 
just the presence of film cooling holes causes transition to 
turbulent flows on the blade surfaces. In later stages, film 
cooling may not be used, and the Reynolds numbers are 
lower. Heat transfer may remain a concern if blades have 
only internal cooling. Hourmouziadis[3] showed nearly an 
order of magnitude decrease in Reynolds numbers between 
the high pressure turbine inlet and the last low pressure 
turbine stage. Predicting the performance of low pressure 
turbines is highly dependent on understanding transition, 
especially at cruise conditions, where laminar separation 
can cause severe loss penalties. 

Surface roughness generally adversely affects blade row 
aerodynamic efficiency. Kind et al.[4], Boynton et al. [5], 
Bammert and Stanstede[6,7] reported decreases in turbine 
efficiencies of up to several points due to surface rough- 
ness. Abuaf et al. [8] and Stabe and Liebert[9] showed 
that just polishing blade surfaces can improve efficiency. 
On the other hand Harbecke et al. [10] showed that pro- 
file loss was not increased until a critical roughness height 
was reached. Boyle and Senyitko[ll] showed that high 
Reynolds number surface roughness doubled vane loss, but 
at low Reynolds numbers roughness improved aerodynamic 
efficiency. Roughness modified the Reynolds number and 
surface location at which separation occurred. To utilize 
surface roughness to improve low pressure turbine aero- 
dynamic efficiency requires understanding of the effects of 
roughness on transition. To predict the efficiency decre- 
ment due to roughness at higher Reynolds numbers also 
requires understanding of the effects of roughness on the 
surface boundary layers. 

There are two approaches to calculating rough surface 
boundary layers. The first is the Discrete Element Method, 
and is ideally suited to situations where the roughness ge- 
ometry is known, and is regular and periodic. This method 
has been used by a number of researchers, and shown to 
give accurate results. Taylor et al. [12] , Hosni et al.[13], 
and Stripf et al. [14] used this method to calculate rough 
surface heat transfer. McClain [15] extended the Discrete 
Element Model to account for random roughness, and Mc- 
Clain et al. [16] presented friction factor and heat transfer 
comparisons with data developed using roughness deter- 


mined from actual blade measurements. The characteris- 
tics of measured turbine blade roughness vary widely, and 
a priori Discrete Element Model calculations need detailed 
information regarding the roughness characteristics. 

The second, and older, approach to calculating rough 
surface boundary layers is to calculate the equivalent 
sand grain roughness. Several authors have proposed 
correlations to obtain the equivalent sand grain rough- 
ness. Among these are Sigal and Danberg[17], Dvorak[18], 
Simpson[19], Dir ling [20], van Rij et al. [21] , and Waigh and 
Kind [22]. Common to all these correlations is the require- 
ment that the roughness image or trace be modeled as 
a roughness geometry, in order to obtain values for pro- 
jected and windward areas. The correlation of Koch and 
Smith [23] only uses a statistical value from the roughness 
trace or image. 

The Koch and Smith [23] correlation provides the sim- 
plest way of determining the equivalent sand grain rough- 
ness from the roughness image. The suitability of this 
correlation can be determined by using the statistical pa- 
rameters of regular roughness geometries. A regular rough- 
ness geometry is one that is suitable for the Discrete El- 
ement Method. It was found that better agreement with 
the experimental data was achieved by a modification to 
the Koch and Smith correlation. 

There are four parts to the work presented herein. The 
first determines a correlation for the equivalent sand grain 
roughness, &s, that is consistent with the experimental 
data of Stripf et al. [24] . The second part determines mod- 
ifications to transition start and relaminarization criteria 
to yield heat transfer predictions consistent with the same 
experimental data. The third part compares heat transfer 
predictions using the proposed modifications with data to 
determine the applicability of the modeling to other data 
sets. The fourth part shows heat transfer and aerodynamic 
loss calculations at a unit Reynolds number consistent with 
current engine operating conditions. Even though both 
the determination of the equivalent sand grain roughness 
and the transition criteria are derived from the same data 
they are not wholly dependent on each other. If the reader 
prefers a different method of calculating /cs, the coefficients 
in the transition criteria could be changed to yield the same 
results with different values for ks- 

DISCUSSION of RESULTS 
Equivalent sand grain roughness 

The roughness geometry used in the test of Stripf et 
al. [24] to measure vane heat transfer coefficients is shown in 
figure 1. In their tests the height of the truncated cones, 
as well as the spacing between cones were varied. The 
dimensions of the roughness geometry is given in Table I. 
For all tests the vane axial and true chords were 53 and 94 
mm. Data were obtained for variations in inlet turbulence 
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Table I. Dimension of roughness elements 


Table II. Roughness characteristics 


Label 

k 

/xm 

y 

/xm 

Ra 

/xm 

Rrms 

/xm 

s k 

ks 

Waigh & 
Kind 
/xm 

Koch & 
Smith 
/xm 

Present 

/xm 

rlOm 

10 

1.8 

2.5 

3.2 

1.65 

48 

15 

37 

rlOs 

10 

1.1 

1.8 

2.7 

2.44 

25 

11 

40 

r20m 

20 

3.5 

5.0 

6.4 

1.65 

96 

31 

73 

r20s 

20 

2.2 

3.5 

5.3 

2.44 

51 

22 

79 

r30m 

30 

5.3 

7.5 

9.6 

1.65 

144 

46 

110 

r40m 

40 

7.1 

10.0 

12.8 

1.65 

192 

62 

147 

r40s 

40 

3.8 

6.3 

10.0 

2.69 

89 

39 

160 

r50s 

48 

4.0 

6.7 

10.8 

2.87 

106 

42 

181 

r80m 

80 

14.1 

20.0 

25.5 

1.65 

384 

124 

293 

r80s 

80 

7.6 

12.7 

20.0 

2.69 

177 

79 

320 


Label 

k 

/xm 

d 

/xm 

V 

/xm 

rlOm 

10 

25 

55 

rlOs 

10 

25 

70 

r20m 

20 

50 

110 

r20s 

20 

50 

140 

r30m 

30 

75 

165 

r40m 

40 

100 

220 

r40s 

40 

100 

300 

r50s 

48 

100 

300 

r80m 

80 

200 

440 

r80s 

80 

200 

600 
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Fig. 1 Arrangement of truncated cone roughness elements. 

intensity and Reynolds numbers. Figure 2 shows their 
measured heat transfer coefficients at an inlet true chord 
Reynolds number of 250,000. As the roughness height in- 
creases, the transition moves forward, towards the leading 
edge, for both the suction and pressure surfaces. Only at 
the highest roughness does the heat transfer appear to be 
fully turbulent over the entire vane surface. An important 
feature of these data is that roughness spacing has a very 
minor effect on heat transfer. 

A number of correlations have been proposed to cal- 
culate the equivalent sand grain roughness, ks- These cor- 
relations show that ks is strongly dependent on roughness 
spacing. These correlations for the ratio of ks to physical 
height are triangular in shape, with a spacing parameter 
that gives a maximum ratio. Closely spaced roughness 
elements are required to achieve the maximum value for 
ks- Stripf et al.[24] gave ks values using the correlation of 
Waigh and Kind [2 2]. This correlation predicted a strong 
effect due to spacing, such that the widely spaced rough- 
ness would have to have nearly twice the roughness height 
to achieve the same value for ks- Figure 2 shows that 
this is inconsistent with the experimental data. The cor- 
relation of Sigal and Danberg[17] has three regions, where 
a plateau separates the two legs of the triangle. At this 
plateau the ratio of ks/k = 8. When applied to rough- 
ness measurements of in service blades, a common feature 
of these correlations is that the roughness image or traces 
have to be transformed into a model geometry. 


r80m 

r40m 

r30m 

r20m 



- 0.08 


- 0.04 0 0.04 

Surface distance, m 


0.08 


0.12 


Fig. 2 Heat transfer coefficients, Re in = 2.5 x 10 5 ,Tu 
8 % 


The correlation of Koch and Smith [23] uses only a sta- 
tistical measure of the surface roughness, without requiring 
that a model geometry be constructed. This correlation 
gives ks = 6.2 R a . Table II gives statistical quantities R a , 
Erms ? and skewness (5k), of the roughness tested by Stripf 
et al.[24]. Also shown are ks values determined from the 
correlations of Waigh and Kind[22], Koch and Smith[23], 
and the present proposed correlation. Since the spacing af- 
fects the value of R a , the correlation of Koch and Smith [23] 
shows a spacing effect not seen in the data. Table II also 
shows that as the spacing increases R a increases, but the 
skewness of the roughness decreases. The three statistical 
quantities given in Table II are defined as: 

r» _^i=i\yi-y\ 

TC a — 

n 

where n equals the number of points in the roughness mea- 
surement, and is generally very large. 
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Fig. 3 Comparison of fully turbulent calculations with 
data. Re m = 250, 000, Tu = 8% 


Rrms = 


and 


£p=ite/i -y ) 2 

n — 1 


siUO/i-y ) 3 

n ^RMS 


s k = 


Since the skewness, S k , is related to Rums, the Koch and 
Smith [23] correlation was modified to give a ks consistent 
with the data shown in figure 2. The revised correlation 
for ks is: 


ks = 4.3#rms( 1 + C'sK-S'k) 


for Sk > 0 


The coefficient of 6.2 was reduced to 4.3 based on the ratio 
of R a to i^RMS- The data in figure 2 and the statistics in 
Table II suggest that Csk is slightly less than one. How- 
ever, for simplicity the constant Csk was set to 1.0. The 
requirement that S k > 0 is not expected to be a difficulty. 
The measurements of Taylor [25] for the roughness of in 
service blades showed positive skewness for many cases. 
Bons et al.[26] measured vane surface roughness, and, ex- 
cept for spalled regions, generally found positive skewness 
values. Spalled regions gave negative skewness. A conser- 
vative approach would be to take ks as the maximum of 
ks = 4.3#rms( 1 + Koch and Smith[23] criterion of 
ks = 6.2 R a , or ks = 4.3Rrms- 

The reasonableness of this correlation is shown in figure 
3 where the data shown in figure 2 are compared with fully 
turbulent calculations for several different equivalent sand 


grain roughness heights using the above expression for ks • 
Since the calculations are for fully turbulent flow, com- 
parisons are meaningful only after transition is complete. 
Overall, the agreement in this region is good. The com- 
parison for the highest roughness, where the measurements 
indicate nearly fully turbulent flow, shows very good agree- 
ment for both the pressure and suction surfaces. 

A question which arises from the data in Table II is 
whether an equivalent roughness height of nearly a third of 
a millimeter is too high compared to actual measurements. 
Most data in the literature gives the roughness height in 
terms of R a , and for comparison with Table II entries the 
Koch and Smith conversion to ks will be used. Tarada and 
Suzuki [2 7] reported a wide variation in R a values. They 
were between 25 and 150/im, ( ks between 155 and 930/im 
), and depended on the engine’s operating environment. 
Taylor [25] measured R & values between 2 and 11/im, (ks 
between 12.4 and 68/im). Bogard et al.[28] extensively 
measured surface roughness on two vanes. The first had 
R & values between 9 and 22/im, and the second ranged 
from 21 to 46/rni. This represents a range of ks between 56 
and 285/rni. Bons et al.[26] measured surface roughness on 
a variety of turbine blades. For deposition and corrosion 
they measured R a values between 3 and 33/im, giving a 
ks range between 19 and 205/im. Zhang and Ligrani[29] 
calculated a ks value of 62/xm for a turbine blade used 
in a utility power application. Just in terms of roughness 
height, the answer is that a roughness height of nearly a 
third of a millimeter is not excessive. 

When Reynolds number effects are accounted for, even 
the highest roughness may not be sufficiently high. Bogard 
et al.[28], in a low Reynolds number experiment, scaled up 
the measured roughness by a factor of 25 to account for 
Reynolds number effects. Since the roughness parameter, 
i7 + , (sometimes referred to as Re k ), increases almost lin- 
early with unit Reynolds number in turbulent flow regime, 
engine unit Reynolds numbers should be matched if the 
roughness height scale factor is one. The data of Stripf et 
al.[24] were obtained at maximum unit Reynolds number 
less than typically seen by the first stage vane at takeoff, 
or in an utility application. On the other hand, the unit 
Reynolds number decreases in later stages, where rough- 
ness transition is more likely to be a factor. 

Heat transfer modeling assumptions 

The data in figure 2 clearly show transition occurring 
for both the suction and pressure surfaces of the vane. 
The locations for transition are strongly dependent on the 
roughness height. Three modeling issues with respect to 
transition were investigated. The first was the start of 
transition. The second was the length of transition. The 
third was relaminarization. A number of approaches were 
examined for each of these modeling issues. Rather than 
detail all of the approaches taken, the modeling that gave 
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the best agreement with all of the data is presented. 

Start of transition. For the start of transition 
Mayle’s[30] transition start model was modified to ac- 
count for roughness effects. In this model when the 
momentum thickness Reynolds number, Re#, exceeds a 
critical Reynolds number, Re grit, transition begins. For 
smooth surfaces Rcqrit = 4007\t -0 - 625 . The critical 
Reynolds number for rough surface transition is: 


Rcrough = 


Rcsmooth 


1 + Tu~ 0 - 625 (c S t{H+ - 5)) 


1.25 


Unless stated otherwise Cst = 0.05. H + is calculated from 
H+ = ksU F sP^/Ci/2/fl The term H+ — 5 occurs because 
a surface roughness that gives H + less than 5 is considered 
hydraulically smooth. Comparisons are made for two in- 
let turbulence intensities, and including the term Tu ~ 0 625 
gave better agreement with data. This term means that 
at high H + values, the start of transition is independent 
of the local turbulence intensity. The roughness term was 
raised to the 1.25 power based on comparisons with the 
data, which show suction surface transition moving from 
near the trailing edge to near the leading edge as the rough- 
ness height increases by a factor of eight. 

The local turbulence intensity used in the transition 
start criteria was calculated using the work of Steelant and 
Dick[31], where the local Tu is calculated from: 


Tu = Tu 1N (U 1N /U FS ) 3/2 

Calculations of the turbulence distribution within the vane 
passage using a k — uj turbulence model was consistent with 
this approximation for the local turbulence intensity when 
the inlet length scale used in the calculations was small, 
Because of the accelerations through the vane passage, the 
maximum local turbulence intensity at the start of transi- 
tion for the smooth vane was less than one percent. The 
data showed that this low turbulence intensity in the tran- 
sition criteria was appropriate. 

Transition length. In the base cases the transition 
length model was not modified to account for roughness 
effects. When modifications were made the spot produc- 
tion parameter, AT, was increased by: 

Trough = Asmooth x exp^CLT(# + — 5)j 

With Cut = 0 there is no modification to the transition 
length. The transition length model is the one detailed by 
Boyle and Simon[32]. It is a modification of the transition 
length model of Solomon et al.[33] to account for Mach 
number effects. It was shown to give good agreement with 
smooth blade heat transfer data for both stator, (Arts et 
al.[l] and Hylton et al. [34]), and rotor test cases, (Arts 
et al. [35] ) . It was found that the length of transition, as 


evidenced by the rough surface heat transfer data, was rea- 
sonably well predicted for many, but not all, cases without 
modifying the spot production parameter, N. 

Relaminarization. Because much of the comparisons 
are at flow conditions where the favorable pressure gradi- 
ents are strong enough to cause relaminarization, it was 
necessary to modify the relaminarization criteria. For 
smooth surfaces it is accepted that when the pressure gra- 
dient parameter, R, exceeds 3 x 10 -6 a turbulent bound- 
ary will relaminarize. Since data show that as the height 
of the roughness increases, the laminar pressure surface 
boundary layer becomes turbulent sooner. Calculations 
without a relaminarization model showed early pressure 
surface transition for many cases. The calculation of the 
critical value of K for relaminarization is: 


( Kcrit)rough 


(Kcrit) smooth 
exp(— C rlm(R + — 5) 1 - 25 ) 


where Crlm = 0.0354, and (Kcrit) smooth = 3.0 x 
10 -6 . The local value for R is calculated from the the 
local inviscid velocity gradient as: 


_ v dUps 
"U 2 S dS 

where S is the streamwise distance from the stagnation 
point. 

Heat transfer calculations. Comparisons are made with 
midspan heat transfer data, so that a two dimensional 
analysis is appropriate. Heat transfer and loss calculations 
were done using a quasi-3d Navier Stokes analysis. The 
solver used was the quasi-3D Navier- Stokes code RVCQ3D, 
run as a two dimensional analysis. This code has been doc- 
umented by Chima [36]. C-type grids, typically 377 x 55, 
were used. Boyle and Simon [3 2] give a more detailed de- 
scription of the analysis. The solutions were monitored to 
assure that convergence was achieved. 

An algebraic turbulence model described by Chima et 
al.[37] was used. It is a two layer model, and incorporates 
the Cebeci- Chang roughness model described in reference 
38. This model increases the mixing length to account for 
roughness. The distance increment is given by: 


A y+ = 0 .9(Vh+ - H+exp-° 167H+ ) 

and 


Ay = Ay + y/(U FS pyfci/2) 

In the algebraic models the increment in y is only applied 
in the inner region. The turbulent eddy viscosity in the 
inner region, // tj j is given by: 


Mt,i = p(dU/dy) 


n{y + Ay)(l — exp (y + + A y + )M + ) 
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Fig. 4 Transition modeling predictions, Re in = 250, 000 


In the two equation k — cu turbulence model roughness 
effects are accounted for by modifying the wall boundary 
condition on cu. Two equation models are most applicable 
when flows are turbulent. Calculations were done using 
the k — uo turbulence model only for fully turbulent cases. 
Calculations using the k — uj model were similar to those 
shown in figure 3, using the algebraic model. For each ks 
value the increase in heat transfer over the smooth calcu- 
lation was nearly the same for both turbulence models. 

When the flows were predicted to not be fully turbu- 
lent, the laminar viscosity was increased to account for the 
effects of high freestream turbulence. Ames et al. [39] de- 
scribed two turbulence models to augment laminar viscos- 
ity. Based on the results presented by Boyle et al.[40], the 
Ames model without a leading edge correction was used. 

Baseline heat transfer comparisons. 

In this section heat transfer comparisons are made 
with the cases used to determine the rough surface heat 
transfer modeling. In a subsequent section comparisons 
are shown for cases that were not used to determine the 
transition modeling. The models for rough surface transi- 
tion and relaminarization have only a few parameters. The 
number is far less than the number of cases that will be 
shown in this section. The goal of the work is reasonable 
agreement for a wide range of cases. Because of the large 
number of cases, and the small number of parameters, it 
is unreasonable to expect perfect agreement for all cases. 

Comparisons of predicted and measured heat trans- 
fer coefficients are given for three Reynolds numbers and 
two inlet turbulence intensities. Figure 4 compares pre- 
dicted and measured heat transfer coefficients for an inlet 
true chord Reynolds number of 250,000 at inlet turbulence 
intensities of 4 and 8 percent. Since the transition mod- 
eling uses only a few parameters, it is note worthy that 
the agreement for the start of transition is good over the 
entire range of roughness heights. The agreement is reason- 
able for both the suction and pressure surfaces, where the 
streamwise pressure gradients are very different. The sec- 
ond highest roughness shows transition close to the leading 
edge for the suction surface, but nearly half way back on 
the pressure surface. 

The analysis under predicts the heat transfer in the 
leading edge region. However, the predictions are in good 
agreement with pressure surface data, where the flow is 
laminar. The under prediction in the leading edge region 
is not due to roughness. Here the analysis under predicts 
the smooth surface heat transfer by more than the small 
increase in heat transfer due to roughness. 

Figure 4 also shows a fully turbulent calculation for a 
smooth vane. The data show that the heat transfer rates 
for the highest roughness can be nearly twice as great as 
those calculated using the smooth fully turbulent assump- 
tion. However, at low to moderate roughness, a fully tur- 
bulent calculation can significantly overestimate the heat 
transfer rate. 
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Comparing figures 4a and 4b shows that the heat 
transfer distributions do not change significantly, even 
when the turbulence level doubles. The critical Reynolds 
number for the start of smooth surface transition decreases 
by nearly fifty percent when the turbulence level doubles. 
By multiplying the roughness term for the start of transi- 
tion by a turbulence intensity factor, the agreement with 
data is improved. However, pressure surface transition at 
the highest roughness is predicted to occur further down- 
stream than the data. This suggests that the exponent on 
the Tu term in the denominator for the start of transition 
should be more negative. 

The question of whether the transition length correla- 
tion should be modified to account for roughness is com- 
plicated by the failure of the transition model to accu- 
rately predict smooth surface transition. For the smooth 
vane suction surface transition occurs in regions where the 
pressure gradient is very negative. At this Reynolds num- 
ber the pressure gradient at the start of transition for the 
smooth surface is more negative than the pressure gradient 
data of Gostelow et al.[41], which was used by Solomon 
et al. [33] to determine the value of the spot production 
parameter, N. The length of transition is very sensitive 
to negative pressure gradients. Figure 4a shows that an 
abrupt smooth suction surface transition was calculated 
for the 8% turbulence level. The transition length would 
increase if N were limited to a smaller value. The data of 
Gostelow and Walker [42] strongly indicate that extrapo- 
lating the correlation of Gostelow et al. [41] to more neg- 
ative pressure gradients would overestimate the value for 
N. At the lower turbulence intensity, figure 4b, the start 
of suction surface transition is at an even more negative 
pressure gradient, because the critical Reynolds number 
for the start of transition is increased. At the lower pres- 
sure gradient the rapid increase in intermittency is appro- 
priate. Fortunately, the question of whether the transition 
model should be modified to account for surface rough- 
ness is not dependent on extrapolating the database for 
the correlation of N. Where rough surface transition is 
seen for negative pressure gradient values, they are within 
the database used to correlate N. The form of the corre- 
lation for N is asymptotic for large positive values of the 
pressure gradient. 

Figure 5 shows the heat transfer comparisons for a 
lower inlet Reynolds number of 140,000 for the same two 
inlet turbulence intensities. Here the analysis shows the 
same trends as the data, and overall is in reasonably good 
agreement with the data. Transition is predicted some- 
what earlier than is seen in the data, especially at the 
higher Tum- Improved agreement with the data would be 
achieved by changing the coefficient in the transition start 
model, Cst, from 0.05 to 0.033. However, the agreement 
with the data in figure 4 for the higher Reynolds number 
would not be as good. 




b) Tu m = 4% 

Transition modeling predictions, Rem = 1.4 x 10 5 
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Table III. Roughness characteristics for Blair’s rotor 
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Fig. 6 Transition modeling predictions, Re in = 0.9 x 

10 5 , Turn = 8 % 


Figure 6 shows the heat transfer comparisons for the 
lowest inlet Reynolds number of 90,000. The comparison 
is at a turbulence intensity of 8%. The agreement is very 
similar to that in figure 5 for the intermediate roughness. 
The results in figures 5 and 6 indicate the 1.25 exponent 
on roughness term should be increased. However, when 
this was done, either in isolation or in combination with 
changing the coefficient, Cst, in the transition start model, 
overall agreement with data was not improved. 

Other rough surface heat transfer comparisons. 

Heat transfer comparisons were made for three other 
data sets, which were not used to deterine the modeling 
for transition. Each data set shows transition behavior in 
the heat transfer results. 


Turbine rotor of Blair [43]. Figure 7 compares Stanton 
number predictions with those measured by Blair [43] at the 
midspan of a turbine rotor. Data was obtained in a large 
scale rotating turbine test facility. Figure 7 shows com- 
parisons at the highest inlet Reynolds number of 580,000. 
Stanton numbers were given for a smooth surface, a near 
smooth surface, and a rough surface. Table III shows the 
roughness measurements from the reference both as ab- 
solute numbers, and as a fraction of axial chord. Only 
Rrms and the maximum-to- minimum roughness heights 
were given. At this Reynolds number the smooth data 
show that the pressure surface is laminar. The predictions 
for the smooth suction surface are lower than the data. 
The inlet turbulence intensity was estimated to be 5%. 


Label 

Rrms 

fim 

^MAX-to-MIN 

fim 

R MAX— to— MIN /C x 
% 

Smooth 

0.33 

7.6 

0.0047 

Near smooth 

6.4 

51 

0.032 

Rough 

N. A. 

660 

0.41 



Fig. 7 Comparison with rotor data of Blair [43] 


Increasing Twin improved agreement for the smooth suc- 
tion surface, but increased pressure surface heat transfer. 

The near smooth data showed transition occurring 
midway along the pressure surface. A ks of 27 gave H + 
values less than 5, so that the predicted heat transfer did 
not change from the smooth surface prediction. Since in- 
creased heat transfer is seen in the data, ks > 4.3Rrms 
should have been used. Perhaps the near smooth rough- 
ness had positive skewness. Interestingly, good agreement 
with the near smooth data is achieved when relaminariza- 
tion is suppressed. 

The rough surface data and predictions show that 
transition occurred close to the leading edge. The rough 
surface ks value was estimated using the near smooth 
Rkms value and the ratio of maximum to minimum 
heights. Clearly, this is a some what crude approxima- 
tion. Nevertheless, the agreement is reasonably good, and 
improves when ks is increased by 50%. A prediction is 
also shown using the k-omega turbulence model. Except 
for the leading edge region, these results agree well with 
the data, and are an improvement for the rear portion of 
the suction surface. 
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Table IV. Rough vane roughness characteristics. 


Trace 

Rrms 
(i m 

R& 

fi m 

Skewness 

ks 

/i m 

1 

17.7 

14.0 

0.206 

92 

2 

20.5 

15.9 

-0.441 

99 

3 

17.5 

13.9 

-0.010 

86 

4 

17.0 

13.0 

-0.085 

81 

Avg. 

18.2 

14.2 


88 


__ C =0.10 C T =0.42 

- <4=0.10 c LT =o 



Fig. 8 Comparison with vane data of Boyle and Seny- 
itko[44] 

Turbine vane of Boyle and Senyitko[44]. Figure 8 com- 
pares measured and predicted Nusselt numbers for the data 
of Boyle and Senyitko[44]. The roughness characteristics 
for four measurements are given in Table IV. Calculations 
were done with ks = 88 fim. The data in figure 8 clearly 
show transition occurring midway along both the pressure 
and suction surfaces. The axial chord was 52 mm. Even 
though the axial chord Reynolds number is over a million, 
and the average Rrms/Cx was 0.00036, the measured inlet 
turbulence was only 1%. Figure 8 shows that at this low 
turbulence intensity smooth surface calculations predicted 
laminar flow for all of the pressure and suction surfaces. A 
smooth fully turbulent calculation gave heat transfer rates 
much greater than seen in the rough surface data over the 
forward portion of both the suction and pressure surfaces. 
The fully turbulent rough surface calculation gave heat 
transfer rates much greater than the data over most of the 
suction and pressure surfaces. Towards the rear of both 


surfaces the agreement between the calculations and data 
is good. This indicates that the ks value determined from 
the roughness statistics is appropriate. 

Rough surface heat transfer predictions using the base- 
line value Cst = 0.05 showed early suction surface transi- 
tion, and no pressure surface transition. Decreasing Cst 
by 20% gave good agreement with the start of suction sur- 
face transition. The variations of ks seen in Table IV are 
about 10%. Decreasing Cst by 20% was the equivalent of 
decreasing ks by 21%. 

Even with Cst = 0.04 the calculated suction surface 
transition length was longer than seen in the data. Increas- 
ing Clt from 0 to 0.071 showed good agreement for the 
suction surface transition length. However, this calcula- 
tion still did not show pressure surface transition consistent 
with data. A Clt = 0.42 was needed to yield pressure sur- 
face heat transfer consistent with the experimental data. 
However, this same calculation gave very poor agreement 
with the suction surface data, due to rapid transition near 
the leading edge. The start of transition was not changed. 
But, suction surface transition began in a region of strong 
favorable pressure gradients. Consequently, the calculated 
transition length with Cur = 0 was very long. 

Turbine blade of Stripf[45]. A second series of tests were 
run with roughness elements similar to those in the tests 
used to determine the roughness transition models. Table 
V gives roughness characteristics for the cases. The cases 
labeled L27a, L27b, and L27c had the same roughness 
heights, but different spacing between the elements. The 
different spacings result in different values of the roughness 
statistics. 

Figure 9 compares the measured and predicted heat 
transfer for this case. Even though these data were not 
used to develop the roughness transition models, the agree- 
ment between the calculations and measured data is good. 
The calculations were done with no modification of the 
transition length due to roughness. The results indicate 
that none was needed. 

Figure 10 compares results for different roughness 
spacing. As seen in Table V the proposed correlation for ks 
has little variation among the three cases with the same 
roughness height. The data in figure 10 shows a depen- 
dency on roughness spacing. As seen in Table V the vari- 
ation of R & values is consistent with the the experimental 
data, in that smaller R & values give lower heat transfer. 
However, just using ks = 6.2 R a does not give good agree- 
ment with the data. The value of H + are too low. 
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Table V. Roughness characteristics of LPT blade 


Label 

k 

fim 

y 

fim 

R a 

fim 

Rrms 

fim 

S k 

k s 

Koch & 
Smith 
fim 

Present 

fim 

L17 

17 

5.2 

6.1 

6.9 

0.83 

38 

77 

L27a 

27 

6.8 

8.6 

10.0 

1.12 

53 

121 

L27b 

27 

3.0 

5.0 

7.5 

2.42 

31 

127 

L27c 

27 

1.7 

3.1 

5.8 

3.54 

19 

131 

L50a 

50 

8.0 

12.5 

16.8 

1.84 

78 

223 

L90 

90 

14.7 

20.8 

27.0 

1.74 

129 

395 
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Surface distance, m 

Fig. 9. Comparison with LPT blade data of Stripf [45] . 

Increasing Cst from 0.05 to 0.10 to adjust the start of 
transition when using the Koch and Smith correlation does 
not give good agreement with the data. It could be argued 
that a coefficient slightly less than one should be applied 
to the skewness term. If this is done the coefficients in 
the roughness models should also be adjusted. Except for 
specialized test cases, roughness typically exhibits notice- 
able variations when measured at different blade surface 
locations. 

Extension to higher Reynolds numbers. 

In typical engine applications the Reynolds number 
may be significantly higher than those used for data com- 
parisons. A series of calculations were done with the same 
physical roughness, but at a higher Reynolds number. 
When the Reynolds number was increased, calculations 
were done for two exit Mach numbers. The lower one, 
(M 2 = 0.55), is the same as for the highest Reynolds num- 
ber tested by Stripf et al. [24] . Since the same geometry 



— kg=395 

— k s = 223 
ks=121 



□ L27a 
O L27b 
L27c 



Fig. 10. Effect of roughness spacing on heat transfer. 



was used, the vane loading moved aft as the Mach number 
increased. The vane isentropic Mach number distributions 
for both exit Mach numbers are shown in figure 11. The 
higher Mach number results in a very aft loaded vane. This 
is a consequence of the vane not being designed for this 
high subsonic Mach number. The results for M 2 = 0.9 
are presented because the distribution is not unlike the 
midspan distribution seen when there is an appreciable de- 
crease in the span between the leading and trailing edges. 
The rapid shock like decrease in Mach number, and con- 
sequently rapid diffusion, is a result of the blade profile. 
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Fig. 12 Nusselt numbers at Re i = 300, 000cm -1 , Tu = 8% 
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Fig. 13 Calculated losses at Re i = 300,000cm -1 

The axial chord Nusselt number distributions for a 
range of roughness heights are shown in figure 12. Nusselt 
numbers are shown for a vane exit unit Reynolds number 
of 300,000cm -1 . The value is less than what is typical 
for a first stage vane of a high pressure turbine. A unit 
Reynolds number in this range is appropriate for the sec- 
ond stage vane, which may or may not be film cooled. The 
unit Reynolds number of 300,000cm -1 is nearly twice as 
great as the highest Reynolds number tested by Stripf et 
al.[24]. For the smooth surface, the distance from the lead- 
ing edge at which suction surface transition is seen nearly 
doubles as the exit Mach number is increased. This is the 
result of the aft loading at the higher exit Mach number. 
At both exit Mach numbers the smooth pressure surface 
remained almost entirely laminar. 

At both exit Mach numbers the effect of increasing sur- 
face roughness is to move the location of the start of suction 
surface transition forward, towards the leading edge. The 
higher exit Mach number, with the aft loading, shows that 
a higher ks value is required for transition to start at the 
same surface location. Except for the smallest roughness 
heights, pressure surface transition is very similar at both 
Mach numbers. Relaminarization is not an issue, and once 
a critical roughness height is reached, transition to turbu- 
lent flow occurs close to the leading edge. For small rough- 
ness heights pressure surface transition behaves similar to 
suction surface transition. 

Because the calculations shown in figure 12 are not 
compared with experimental data, it is useful to relate 
these results to those presented in the literature in a qual- 
itative manner. This is done to show the reasonableness 
of the modeling presented herein. Abuaf et al.[8] exam- 
ined the effects of small roughness on vane heat transfer 
and aerodynamic efficiency. Their exit vane axial chord 
Reynolds number varied between 6.8 x 10 5 and 2.3 x 10 6 . 
The inlet turbulence intensity was 14%. Their maximum 
R & value was 2.33/zm, but the axial chord was only 4.8cm. 
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The vane was forward loaded, and at the lowest Reynolds 
number suction surface transition occurred just upstream 
of the peak isentropic Mach number. Roughness was not 
a factor in in transition at the lowest Reynolds num- 
ber. As the Reynolds number increased, the vane with 
R a = 2.33/im, ( ks ~ 15), was seen to move transition 
closer to the leading edge. At the highest Reynolds number 
transition was seen at the leading edge. Pressure surface 
transition was not affected by roughness, but was turbulent 
for nearly 80% of the surface distance. In the turbulent re- 
gion the heat transfer rates increased approximately 10 to 
15 percent as a result of surface roughness. 

The importance of roughness in terms of aerodynamic 
losses is illustrated in figure 13. The loss coefficient, Y is 
calculated from: 

Y _ ~ 

Pin, t — P2 

Figure 13 shows the loss coefficient, F, as a function 
of ks- At the highest equivalent height, the loss coefficient 
nearly doubles, and approaches an asymptotic value. The 
increase in smooth surface loss for the higher exit Mach 
number remains even as the surface roughness increases. 

The data show that an important parameter is H + . 
In fully turbulent flow H + is nearly proportional to the 
Reynolds number. Therefore, doubling the unit Reynolds 
number results in a similar loss distributions, but with the 
abscissa in figure 13 reduced by half. 

CONCLUDING REMARKS 

An approach to calculating an equivalent sand grain 
roughness has been given. This approach is consistent with 
experimental data in that the weak effect of roughness 
spacing is accounted for. This approach uses only statisti- 
cal quantities, which are readily available from roughness 
measurements. When the calculated ks values were used 
with either an algebraic or a k — u turbulence model, heat 
transfer rates were in good agreement with data in the 
fully turbulent regions. 

Two-dimensional calculations were done to determine 
criteria for transition and relaminarization that were rea- 
sonably consistent with experimental data. Both calcula- 
tions and data show heat transfer rates much less than for 
fully turbulent flow at small roughness heights. For high 
roughness heights heat transfer rates were nearly double 
those for the fully turbulent smooth flow. This simplified 
approach gave reasonable agreement with data for most 
of the cases examined. There were two prominent excep- 
tions. Both involved the blade pressure surface. The data 
for the near smooth case of Blair showed that a roughness 
with a calculated H + < 5 appeared to cause a laminar 
boundary to become turbulent. In the second case a vane 
tested with low inlet turbulence exhibited rapid transition 


midway along the pressure surface. This was not predicted 
by the analysis, and further work in this area is needed. 

Using the modification to account for roughness on 
the start of transition, calculations were then done for a 
unit Reynolds number of 300,000cm -1 . This is represen- 
tative of the second stage of a high pressure turbine. For 
good cycle efficiency the second stage vane might not be 
film cooled, and accurate transition predictions would be 
important. Vane loading was seen to be an important pa- 
rameter. Aft loading delayed transition, both for smooth 
and slightly rough surfaces. Even at this moderately high 
Reynolds number, smooth and slightly rough surfaces had 
heat transfer rates significantly less than for the fully tur- 
bulent calculations. On the other hand, high roughness 
levels showed heat transfer rates far greater than those 
calculated for fully turbulent flow. 

The effects of surface roughness on losses were also 
examined. The smooth aft loaded vane had a greater 
predicted loss than the forward loaded vane. Loss lev- 
els approached asymptotic values nearly twice that of the 
smooth vane as ks increased. 
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